Reformation of an oblique shock observed by ClusteiQ 
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On 16 March 2005, the Cluster spacecraft crossed a shock almost at the transition between the 
quasi-perpendicular and quasi-parallel regimes (Obti ~ 46°) preceded by an upstream low- frequency 
(~ 0.02 Hz in the spacecraft frame) wavetrain observed for more than 10 mn. The wave semi- 
cycle nearest to the shock was found to grow in time, steepen and reflect an increasing fraction 
of the incoming ions. This gives strong indication that this pulsation is becoming a new shock 
front, standing ~ 5Ap upstream of the main front and growing to shock- like amplitude on a time- 
scale of ~ 35Qp. Downstream of the main shock transition, remnants of an older front are found 
indicating that the reformation is cyclic. This provides a unique example where the dynamics 
of shock reformation can be sequentially followed. The process shares many characteristics with 
simulations of reforming quasi-parallel shocks. 



I. INTRODUCTION 

The question of the stability of the shock structure 
has regained recent interest with the advent of mul- 
tipoint measurements at the Earth's bow shock. Nu- 
merical simulations have indicated for a long time that 
even under steady upstream conditions and for a broad 
range of parameters coUisionless shocks can exhibit 
dramatic structural changes and eventually self-reform 
quasi-periodically [2l| . However experimental ev- 

idence of this kind of behaviour in space plasmas have 
been limited. For instance, Thomsen et al. [38| and 
Thomsen et al. [39| have observed quasi-periodic varia- 
tions on a time-scale of approximately two upstream pro- 
ton gyroperiods in the ion velocity distributions down- 
stream of the quasi-parallel Earth's bow shock which 
they have attributed to shock reformation. Using multi- 
point Cluster measurements, Lucek et al. have stud- 
ied the scales and growth of large amplitude pulsations 
(Short Large Amplitude Magnetic Structures, SLAMS) 
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thought to be part of the quasi-parallel shock transition 
region. At the quasi-perpendicular shock, Horbury et al. 
|la | have noticed on some Cluster crossings a significant 
variability in magnetic field profiles, particularly in the 
foot, despite a relatively small spacecraft separation. A 
large variability in magnetic and electric field time se- 
ries was also observed by Lobzin et al. [2^ at a high 
Mach number quasi-perpendicular shock crossing, com- 
plemented with bursty variations in reflected ions occur- 
ing on a time scale comparable to a proton gyroperiod. 
Another form of non-stationarity was found by Moullard 
et al. [2^ who have identified coherent oscillations with a 
wavelength of a few tens of upstream ion incrtial length 
confined to the shock layer and propagating along it. 

The general problem of the shock stationarity is com- 
plicated by the importance of the shock geometry. Quasi- 
parallel and quasi-perpendicular shocks are indeed dis- 
tinguished by significantly different structure, scales and 
dissipative processes (e.g. [11]). At least for supercrit- 
ical shocks resistivity is insufficient to provide the re- 
quired dissipation, which is then provided in a first step 
by refiecting upstream part of incoming ions. Assum- 
ing a specular refiection process and a step-like shock 
transition one finds out that the limit between quasi- 
perpendicular and quasi-parallel shocks Obu ~ 45° cor- 
responds to the angle below which the guiding-centre of 
the reflected ions is directed upstream instead of being di- 
rected back towards the shock [13]. However specularly 
reflected ions at shocks only slightly below Obu = 45° 
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may still re-encounter the shock front during the course 
of their gyromotion unless the shock angle is actually 
lower than ~ 40° [3l| . Only then are these reflected ions 
able to escape upstream, and by backstreaming against 
the solar wind flow they excite upstream waves (e.g. 
[131 )■ As these low- frequency waves grow to large ampli- 
tudes they interact with the incoming solar wind, giving 
quasi-parallel shocks a more extended and visibly more 
complex and dynamic structure than quasi-perpendicular 
ones. Clearly there is more to coUisionless shock dissipa- 
tion than simple specular reflection. In the case of curved 
shocks such as the Earth's bow shock this picture is fur- 
ther complicated by ions streaming from other parts of 
the shock which also contribute to instabilities near the 
quasi-parallel shock. Burgess et al. fi\ provide a recent 
review of observations of quasi-parallel shock structure 
and processes (see also [^). 

Observations have therefore shown that the bow shock 
can significantly deviate from the textbook picture of a 
locally planar and stationary structure and that it may 
eventually reform. However the scales and complexity 
of shock reformation combined with the limited number 
of measurement points have not provided so far direct 
evidence by following sequentially the reformation of a 
shock front. 

This paper presents a case study of a crossing by the 
Cluster fleet of a shock with Osn ~ 45°. The region 
upstream of the shock exhibits a low-frequency quasi- 
monochromatic wave as more commonly observed up- 
stream of quasi-parallel shocks. The wave cycle nearest 
to the shock ramp is found to steepen and grow into a 
pulse-like structure which we argue corresponds to the 
formation of a new shock-ramp. The four Cluster space- 
craft were able to sequentially observe this process which 
we interpret to be shock reformation. 



II. DATA 

This study concentrates on a shock crossing on 16 
March 2005 around 1530 UT by Cluster occurring dur- 
ing a data burst mode interval. The quartet was in a 
tetrahedron configuration with inter-spacecraft separa- 
tions of the order of 1300 km. The magnetic field data 
from the Flux-Gate Magnetometer (FGM) Q used in 
this study was averaged at a resolution of one vector per 
spin-period (4 s) or 10 vectors per second, averaged from 
a 67 vectors/second sampling frequency in burst mode. 
The electric field data from the Electric Field and Wave 
(EFW) JJj instrument used here has a temporal resolu- 
tion of approximately 2 ms. Ion data was provided by the 
Hot Ion Analyzer (HIA) of the Cluster Ion Spectrometer 
(CIS) [28J which measures fluxes of positive ions irre- 
spective of species in the energy range 0.005-26 keV/e 
and takes a spin period to build a full 3d distribution 
(transmitted to the ground every spin in burst mode) 
and has an angular resolution of 22.5° x 22.5°. How- 
ever during this interval the instrument was in a mag- 



netospheric mode not optimized for the solar wind, and 
CIS/HIA data were available on spacecraft 1 and 3 only. 
Electron data came from the Low Energy Electrostatic 
Analyzer (LEEA) of the Plasma Electrons And Current 
Experiment (PEACE) [17(]. In this burst mode, a 3d 
distribution with reduced polar resolution (3DX1) was 
transmitted to the ground every spin, and consists of 
26 energy levels in the range 0.007-1.7 keV, 6 polar and 
32 azimuthal angular bins. This data was then reduced 
on the ground to pitch-angle distributions using high- 
resolution FGM data, and corrected for the spacecraft 
potential effect using spin-resolution EFW data. Mag- 
netic field and particle from the MAG and SWEPAM in- 
struments onboard the Advanced Composition Explorer 
(ACE) spacecraft fs^l were used to compute the upstream 
solar wind parameters. The time delay taken by the so- 
lar wind to travel from the Lagrange point to the Earth's 
bow shock was taken into account assuming it moves with 
a constant velocity parallel to the GSE x-axis. 



III. SHOCK OBSERVATIONS 

A. Global shock properties and Cluster 
configuration 

The solar wind as monitored by ACE, time-shifted by 
the solar wind travel time from the Lagrange point to the 
shock, remains steady and quiet during the whole time 
interval and for nearly half an hour before, in particu- 
lar with a velocity around 400 km/s and no important 
change in magnetic field direction. Because of the pres- 
ence of large amplitude fluctuations extending far up- 
stream of the shock (and of the CIS instrument mode), 
the asymptotic upstream field and particle parameters 
are estimated from the ACE measurements (taking into 
account the solar wind travel time to the shock). 

The crossing occurs at (12.7, 0, 4.6)i?e (GSE), from 
downstream to upstream (outbound crossing). The 
shock timing analysis (which assume a planar surface 
moving at constant speed) (32j yields a normal n = 
(0.93, -0.12, 0.35) (GSE), a shock angle between the up- 
stream magnetic field and the normal Obu = 45° and 
a velocity along normal of -13 km/s in the spacecraft 
frame. The Abraham-Shauner method (32j which makes 
use of the MHD jump (Rankine-Hugoniot) conditions for 
the magnetic and velocity fields (and assumes as well a 
planar surface and shock stationarity) yields a similar re- 
suh, n = (0.93, -0.10,0.35) and Obu = 46.5°. The shock 
is therefore within experimental errors at the formal limit 
of quasi-perpendicular and quasi-parallel shocks, and we 
shall generically call it an oblique shock. 

The projected spacecraft locations onto the coplanarity 
plane show that Cluster- 1 (CI), 3 and 4 are approxi- 
mately aligned along the shock normal while CI and 2 
are perpendicular to it and cross the ramp nearly simul- 
taneously (fig.[I]). The shock is crossed first by C4, then 
by CI and 2 together and finally by C3. 
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TABLE I: Main shock parameters at the asymptotic upstream and downstream locations. Upstream parameters were taken 
from ACE data. Frame-dependent quantities are given in the Normal Incidence frame (the shock rest frame where the upstream 
solar wind velocity is directed along the shock normal). 

Parameters units upstream downstream 
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Convected proton gyroradius, Upi/Slgp 
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47 
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Alfven velocity, ca 
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77 
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Alfven Mach number. Ma 




5.5 


1.00 


Magnetosonic Mach number, Mma 




4.6 


0.8 



The shock is super-critical with an Alfven Mach num- 
ber Ma = 5.5. The proton thermal to magnetic pressure 
ratio for this shock is low, /?p = 0.1. Other shock and 
plasma parameters are summarized in Table HI 



B. The upstream low- frequency wave 

An upstream ultra-low frequency wave (ULF) with a 
period ~42 s (/ = 24 mHz) in the spacecraft frame 
is observed for over 10 minutes from the ramp (fig. 
(2). The time-series appear quasi- monochromatic at 4 s- 
resolution, but considerable broadband higher-frequency 
fluctuations are seen in the higher-resolution data. The 
ULF oscillations can be seen in the magnetic field inten- 
sity as well as all three components and strongly increase 
in amplitude in the vicinity of the ramp. For instance, 
at the shock ramp the magnetic field fiuctuations along 
the shock normal reach 5Bm /Bm w 1.7, and the fiuctua- 
tions in the non-coplanar component of the magnetic field 
5Bm are nearly as large as the overall difference between 
upstream and downstream magnetic field component in 
the maximum variance direction B^. The general shock 
profile, ULF wave characteristics and its relative loca- 
tion with respect to the shock seem remarkably similar 
on all four spacecraft despite the large spacecraft separa- 
tion and temporal spread of the crossings (the first and 
last shock crossings are approximately 80 s apart). This 
gives to the wave an unexpected appearance of phase- 
stationarity with respect to the shock (to be discussed 
in the final section). Oscillations extend over one or two 



periods downstream of the ramp with a shorter period 
~36 s, although it is not entirely clear whether they 
correspond to a transmitted wave or to the overshoot- 
undershoot cycle. 

Minimum Variance Analysis (MVA, see e.g. fs^) ap- 
plied suggests that the wavevector is roughly aligned 
with the magnetic field (the deviation, 9kB ~ 24°, is 
mainly in the out of coplanarity plane direction), k//c = 
±(0.63,0.38,0.67) (NML). The wave has a high degree 
of polarization (0.92), is left-hand polarized (with re- 
spect to the magnetic field) in the spacecraft frame with 
ellipticity e = —0.81 and has a low compression ratio 
Ce = [Sn/nof (Bi/SB^) = 0.12 (tab. 

More properties can be derived from the time delays 
between spacecraft assuming a planar wavefront. Al- 
though in general for a monochromatic wave the time 
delays can only be determined modulo the wave period, 
we find that the smallest delays determined from the 
time-series provide a direction close to that of the MVA, 
namely k/k = (0.66, 0.37, 0.65) in NML coordinates (2° 
away from the MVA result). The corresponding phase 
velocity in the spacecraft frame is 142 km/s which yields 
a wavelength of 5930 km. From the wavevector and 
the mean upstream plasma velocity, the frequency in the 
plasma frame is found to be 0.01 Hz and the plasma frame 
phase velocity is 56 km/s, lower than but comparable to 
the upstream Alfven speed (76 km/s). In this frame, the 
wave propagates against the solar wind flow and there- 
fore is right-hand polarized. The properties summarized 
in table |lT] such as wavelength Re, upstream prop- 
agation direction and right-hand polarization are fairly 
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FIG. 1; Shock crossing configuration. Top panels show the 
upstream magnetic field lines and the location of the crossing 
(black square) in the GSE xy and xz planes. The main pan- 
els show the asymptotic magnetic fields, solar wind velocity 
in the normal incidence shock frame, and projection onto the 
coplanarity (NL) plane of spacecraft locations and direction 
of the upstream wavector. The grey area corresponds to the 
main shock ramp, and distances are normalized to the up- 
stream proton inertial length \p. Spacecraft in this plane are 
approximately aligned along the shock normal for CI (black), 
3 (green) and 4 (blue), while C2 (red) crosses the shock nearly 
simultaneously to CI but approximately 10 Ap away along the 
shock front. The spacecraft travel from downstream to up- 
stream. 

typical of ULF waves studied using ISEE [ll] or Cluster 
[l|, II] datasets, and are thought to result from an elec- 
tromagnetic ion/ion right-hand resonant instability due 
to ions backstreaming from the shock. Finally, the phase 
velocity along the shock normal in any shock rest frame 
is estimated to be v^^'^°'^'^-* • fi « —80 km/s showing that 
as expected the wave is not phase-standing but convected 
towards the shock by the solar wind. 

Applying MVA to shorter time intervals (2 wave peri- 
ods) however reveals that these properties change closer 
to the shock front. As shown in fig. |3] the wavevector 
approximately aligns itself with the shock normal, the 
polarization becomes less circular and the compression 
ratio increases. These changes appear during intervals 
which do not yet include the shock ramp, corresponding 
to specific properties of the cycle nearest to the ramp 
which may be affected by the shock foot and reflected 
ions. Since the wave is assumed to be planar and has a 
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FIG. 2: 4s-resolution magnetic field intensity (top panel) and 
components in shock normal coordinates on all four space- 
craft. Data are time-shifted to allow comparisons of the shock 
structure and upstream wavetrain measured by the space- 
craft. Time information is translated into distance assum- 
ing the whole structure travels at the constant shock velocity 
estimated from the timings. 

high degree of polarization, an effect of the alignment of 
the wavevector with the shock normal should be to di- 
minish the perturbations due to the wave to the planarity 
of the shock surface. In addition the oscillations seem to 
remain in-phase with the shock as if the ramp was part 
of one cycle and other cycles were phase- standing next to 
it. Indeed, the timings of the pulse nearest to the shock 
ramp confirm the wavevector alignment with the normal 
(consequently O^b ~ 40°) and a velocity nearly identical 
to that of the shock, about -13 km/s in the spacecraft 
frame. This yields a wavelength an order of magnitude 
lower than further upstream, A ~ SOOfcm « 5Ap. 

Besides changes in wavevector and velocity, the wave 
experiences a strong amplification near the shock. The 
cycle standing nearest to the ramp indeed nearly reaches 
shock-like amplitudes and displays an interesting be- 
haviour detailed in the next section. 



C. The growing and steepening pulse upstream of 
the shock ramp 

As noted in the previous paragraph, upstream mag- 
netic field fluctuations reach their largest amplitude at 
the wave cycle nearest to the shock front. This large am- 
plitude pulse-like structure is crossed by the spacecraft 
in the same order as the shock ramp (C4 first, followed 
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TABLE II: Properties of the upstream ultra-low frequency wave in spacecraft and plasma frames. 

Spacecraft Plasma 
Frequency / fgp 0.14 0.07 
Wavelength /Ap 62 62 
Phase velocity /ca 1.8 0.74 
Polarization degree 0.92 0.92 
Polarization Left-hand Right-hand 

EUipticity -0.81 + 0.81 

Compression ratio 0.12 0.12 




15:28 15:30 15:32 15:34 

i6Mar05 Universal Time [h:m] 



FIG. 3: Wave properties in spacecraft frame derived from 
MVA on 84s intervals of 4s resolution data. The top panel 
shows the magnetic field intensity. Next panel shows the polar 
angle 6 between k and n, and which is the azimuthal angle 
in the LM plane with respect to the L axis. The dotted lign 
indicates Obu- Lower panels show the ellipticity and electron 
compression ratio. 



nearly simultaneously by CI and C2 and then C3), show- 
ing that it is not a partial recrossing of the shock front 
but a distinct upstream structure. The time delay be- 
tween the crossings of the shock ramp and the feature is 
nearly the same for all four spacecraft (about 20s), sug- 
gesting that this structure extends parallel to the shock 
ramp and remains at a constant distance from it. 

As shown in fig. |4] the magnetic field intensity of the 
structure is the lowest at the crossing by C4, larger at 
CI and C2 (and slightly more so at C2 than CI) and 
largest at C3 where its amplitude is comparable to that 
of shock itself. The same observation applies to Bl ■ The 
structure is therefore growing in time. An exponential 
fit to the peak amplitudes yields a growth rate of 7 ~ 
4-10~^ Hz w 0.03/gp„. The amplitudes are equally well fit 
by a linear curve with slope 0.011B„ s^^ « O.OSBufgpu- 
Both models estimate that it takes up to ~ 35 upstream 



— ci 

-1 — C2 

Slope: 0.0042 s - 0.03 / — C3 




15:29 15:30 15:31 

i6Mar05 Universal Time [li:m] 

FIG. 4: Main shock ramp and upstream pulse, shown with 4 s 
resolution FGM data. The pulse growth rate is approximately 
0.03/cp 



proton gyroperiods for the pulse to grow to shock-like 
amplitudes. 

High-resolution magnetic field data show that besides 
growing in amplitude the structure is found to steepen 
(fig. [5]). The steepening is most visible in Bl and occurs 
on the upstream edge of the pulse (apart perhaps from 
C2 on which the "downstream" side of the pulse seems 
at least as steep as its upstream side). As the steepening 
proceeds quasi-periodic whistler-like fluctuations (at « 
0.15 Hz, best seen on Bm) are emitted upstream and the 
measured electric field magnitude (E^ + EyY^"^ increases 
too and reaches on C3 values comparable to those in the 
ramp. Furthermore, high-resolution data reveal that the 
pulse is seen slightly earlier on C2, consistent with its 
position slightly upstream of CI as sketched in Fig 1. 
However, that data also shows the pulse's growth and 
steepening to be more advanced at C2. This shows that 
at this separation scale the pulse's structure and growth 
is not perfectly homogeneous along the shock plane. 

Finally one notes on CI and C2 in between the pulse 
and the shock ramp localized dips in B^ reaching nega- 
tive values, similar to that observed during the reforma- 
tion cycle in Id 0, iiOj and 2d (s^j quasi-parallel shock 
hybrid simulations. 
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FIG. 5: High-resolution (10 vectors per-second) magnetic 
field measurements of the growing structure in shock NML 
coordinates. The structure is not only found to grow in 
amplitude but also to steepen while emitting whistlers, and 
measured electric field intensities correspondingly increase to 
reach ramp-like values. 



D. Particles and cross-shock potential 

If the growing feature is indeed becoming a new shock 
front, then it should affect the heating and reflection of 
incoming solar wind particles and correspondingly de- 
velop an electric potential jump. fig. [6] shows magnetic 
intensity profile and ion energy spectra (from CIS/HIA, 
all species and summed over all view angles) on CI and 
C3. Two distinct ion populations are clearly observed in 
the solar wind. The lowest energy one (« 1 keV) is the 
incoming solar wind beam which undulates under the ef- 
fect of the wave, yielding large oscillations of the velocity 
moment. It seems however that the wave results in little 
or no ion heating (which is nevertheless difficult to check 
quantitatively in the absence of reliable temperature es- 
timates due to the particular CIS instrument mode). 

There is also a higher energy (« 3 keV), less dense but 
hotter population which is strongly modulated by the 
wave suggesting that some of these energetic ions could 
be trapped by the large-amplitude wave. The highest 
count rate of energetic ions on CI appears during the 
shock ramp crossing around 15:28:40 UT. This distinct 
ion population has an energy around 3 keV and should 
consist of gyrating specularly-reflected protons. A simi- 
lar but lower count-rate group of ions is observed on the 
next peak of magnetic intensity. The situation on C3 
is slightly different however. There is no distinct ener- 
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FIG. 6: Energy spectrograms of ions from CIS/HIA, summed 
over all angles. Upstream, both solar wind and energetic ions 
are strongly affected by the wave. At the shock, a distinct 
population of reflected ions around 3 keV is observed at the 
shock ramp on CI but not on C3, where it is observed on the 
upstream steepened structure instead. 



getic population of energetic gyrating ions observed at 
the "old" shock ramp. These are only seen at the up- 
stream steepened structure, as if it became a new shock 
front where most of the ion reflection occurs. One notes 
however that there is little appreciable ion heating be- 
tween this structure and the old ramp, suggesting that 
some solar wind plasma is not processed by a full ramp 
structure but caught in between. 

Specular ion reflection being associated to the cross- 
shock potential at least part of the potential jump must 
occur at the new front. One of the most reliable ways 
to estimate the electric potential across a shock is to use 
electron data combined with Liouville mapping ([2^ and 
references therein). Based on assumptions of conserva- 
tion of electron energy in the de Hoffmann- Teller frame 
and flrst adiabatic invariant, this technique relates the 
changes in electron velocity distributions between two lo- 
cations to the unknown potential difference. A variant of 
this technique is used which takes into account the held 
maxima in between the two locations ^2Q] . 

The estimated cross-shock potentials on all four- 
spacecraft in the de Hoffmann- Teller frame are shown 
in flg. [71 where time-series have been shifted in order to 
make the main magnetic ramps coincide. Their down- 
stream values are remarkably similar on all four space- 
craft despite the shock non-stationarity. However the 
location of the main potential jump with respect to the 
main magnetic ramp seems to differ from spacecraft to 
spacecraft. During the C4 crossing, when the new ramp 
was just starting to form, the potential jump occurs at 



7 




-20 -10 10 20 
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FIG. 7: Magnetic field intensity (top) and estimated electric 
potential (bottom) time-shifted in order to make the magnetic 
ramps coincide. The downstream cross-shock potential values 
are very similar on all four spacecraft. However the location 
of main potential jump seems to move upstream of the main 
magnetic ramp. 



the magnetic ramp. It is then observed further upstream 
on CI and C2 crossings, and near the newly formed ramp 
on C3. 



E. Large-amplitude downstream perturbations and 
indications of a reformation cycle 

Large perturbations are also found downstream of the 
main shock ramp as shown in fig. [8l localised at ap- 
proximately the same distance from the main ramp as 
the upstream pulse. On all four spacecraft a depression 
in the mean magnetic field intensity is observed. Large 
magnetic field and electric field fluctuations are present 
within this depression. The amplitudes of the electric 
field fiuctuations are comparable to or even larger than at 
the main ramp. The depth of the depression and ampli- 
tude of the fluctuations are also decreasing in time, being 
the lowest on C3 which is the last to cross the shock. The 
decrease occurs on time scales comparable to the growth 
of the upstream pulse. These perturbations can be inter- 
preted as remnants of an older shock ramp which decays 
in time, showing that the previously described formation 
of a new ramp is not an isolated event but part of a quasi- 
periodic reformation cycle. Based on the growth rate of 
the upstream pulse and the decay of downstream per- 
turbations, the reformation period can be estimated as 
several tens of upstream proton gyroperiods, and could 
be a few periods of the upstream low-frequency wave. 
This is signiflcantly longer than the period of variations 
of the downstream ion populations "^fpu) observed 
by Thomsen et al. (ssj , although the shocks studied by 



these authors have slightly different parameters and in 
particular higher Mach number. 



IV. SUMMARY AND DISCUSSION 

The shock analysed in this paper is oblique {Obu ~ 
45°), moderate Mach number {Ma ~ 5.5) and low-/? 

(A « 0.1). 

Upstream of the shock a long-wavelength (A w 
62Ap), low- frequency (/ « 0.07/p in the plasma frame) 
and right-hand polarized in the plasma frame quasi- 
monochromatic wave propagates against the solar wind 
flow approximately parallel to the upstream magnetic 
held [OkB ~ 24°). Its properties are consistent with a 
magnetosonic-like wave excited by a weak (rif, <C n) and 
cool {vTb < w&) ion beam backstreaming from the shock 
through an electromagnetic ion/ion right-hand resonant 
instability [TTj . Such an instability has its maximum 
growth rate for field-aligned wavevectors and obeys the 
linear dispersion relation in the cold plasma approxima- 
tion uj = k\\Vb — fip (real part). In the low frequency 
limit Lo VLp one therefore has fc|| ~ flp/vf,. Assum- 
ing that the energetic ion population seen in fig. [5] (^^3 
keV in the spacecraft frame) corresponds to the back- 
streaming ions yields an estimate for the phase velocity 
w/fc|l sa {uj/flp)vf) « 80 km/s in good agreement with 
the timing results, but the wavelength is found to be 
A « 8300 km which is larger than the observed wave- 
length of 5900 km. It is therefore likely that the 3 keV 
population (which actually spreads from 2 to 10 keV) is 
modulated by the wave but not its source. Instead the 
wave may have been generated further upstream by a 
beam of ions coming from another part of the shock. In- 
deed the wave properties described in this paragraph are 
also similar to those of ULF waves in the ion foreshock 
[1) waves which are generally not observed simul- 
taneously to the field-aligned beams but in association 
to diffuse-ions or g yrating beams in cyclotron resonance 
with the wave |24| . The observed high energy ion popu- 
lation modulated by the wave is therefore more likely to 
be either a gyrating beam of refiected ions which man- 
aged to escape from the shock or gyrating ions trapped 
by the wave. 

Due to its small phase velocity the wave is swept back 
by the solar wind, and closer to the shock the wave prop- 
erties change significantly. The wavevector aligns itself 
with the shock normal, which could be due to an increase 
of refraction index in the foot region. A similar align- 
ment was found in the 2D hybrid simulations of Scholer 
et al. [30j . which they attributed to an increase in den- 
sity of the diffuse ions near the shock. This alignment 
makes the perturbations of the shock due to the wave 
more homogeneous along the shock surface. It is this 
greater deviation from parallel propagation which allows 
the wave to become more compressional and steepen, an 
effect which was conjectured by Hada et al. [14] to happen 
as ULF waves are convected into areas of the foreshock 
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FIG. 8: Time-shifted high-resolution magnetic field intensity and electric field (|-Ba;t/| = {E^ + Ey)^^^). Spacecraft are ordered 
in the order by which they cross the shock. Besides the growth of the upstream pulse which becomes a new ramp, data show 
strong perturbations approximately one wavelength downstream of the main ramp. This can be interpreted as remnants of an 
old shock front from a previous reformation cycle. 
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FIG. 9: Illustration of the wave refraction near the shock. 
The wavefronts (dashed lines) become nearly parallel to the 
shock surface (thick grey line), which requires a lower phase 
velocity in the shock frame at location 1 (near the shock) than 
at 2 (further upstream). 

with higher densities of diffuse ions (see also UM)- Be- 
sides, the wave has all the appearance of phase-standing 
with respect to the shock. From a velocity identical to 
that of the shock in this region, the wavelength is found to 
be 10 times shorter than further upstream. This slowing 
down of the wave in the shock frame is entirely consistent 
with its wavevector alignment with the shock normal, as 
illustrated by fig. [9] 

In addition to refractive effects an instability could be 
operating near the shock front, partly accounting for the 
changes in wave properties and its strong amplification. 
This could involve, as dicussed by Winske et al. pio| . 
a relatively dense and cold beam of specularly reflected 
ions observed near the shock front during the reforma- 
tion process in the simulations, or an instability localized 
at the shock front resulting from the interaction between 
upstream and some of downstream heated ions. This " in- 
terface instability" generates waves on the magnetosonic 
branch with wavelength intermediate between whistlers 
and ULF waves. For Ma = 5 the maximum growth rate 
corresponds to kXp « 0.5 or A sa AnXp. This wavelength 
(A « 5Ap) estimated near the shock is smaller, but still in 
the range where the growth rate is close to its maximum 
value. The instability is weakly dependant on Obu and 
tends to have its wavector aligned with the shock nor- 
mal. However the exact mechanism(s) operating on the 
wave near the shock front is (are) certainly more com- 
plicated since we are dealing with very large-amplitude 
(nonlinear) waves and a highly inhomogeneous situation. 
One can however note that because the wavelength near 
the shock is about 6 times the theoretical wavelength 
of the phase-standing whistler the mechanism leading to 
reformation is distinct from the whistler-mediated ref- 
ormation discussed for oblique and quasi-perpendicular 
shocks by Biskamp and Welter Krasnoselskikh et al. 
p^ . Scholer and Burgess [29| . 

Due to the strong amplification of the wave near the 
shock, its last cycle before the shock acquires an am- 



plitude comparable to that of the shock itself. This 
pulsation grows in time, steepens and in the process it 
emits whistlers, reflects more and more ions and con- 
tributes more significantly to the cross-shock potential 
jump. Reasons for this growth and steepening can be the 
change in refractive index, a localized instability near the 
shock front such as the interface instability or nonlinear 
dynamics of the pulse and its interaction with the shock. 
Furthermore the amplification mechanism may not oper- 
ate steadily in time. The negative dips in Bj^ seen only on 
CI and C2 during the intermediate growth stage of the 
pulse may indeed be the signature of a transient dense 
ion beam reflected from the shock at this stage of the 
shock reformation (e.g. as in the simulations of Thomas 
et al. [131). At least this indicates an unusual and tran- 
sient processing of ions in between the main ramp and 
the pulse during the pulse growth. 

Whatever mechanism operates here, these observations 
give strong evidence that the wave cycle standing near- 
est to the shock is forming a new shock front on a time 
scale ^ or less, comparable to that observed in hy- 

brid simulations of quasi-parallel shocks (e.g. [5l. Isol. [37| ) . 
The reformation scale (distance between the shock front 
and pulsation) is 5Ap, also in good agreement with 
simulation results. Assuming that the real shock front 
thus jumps forward such a distance in a time of 
its velocity is found to be only a few km/s larger than 
that of the "former" shock front in the solar wind frame, 
which does not induces a significant change in effective 
Mach number. At approximately the same distance from 
the main shock ramp but on the downstream side there 
are magnetic and electric field fluctuations with ampli- 
tudes comparable to those at the main ramp. These can 
be naturally interpreted as remnants of an older ramp 
which decays on a time scale of at least a few wave pe- 
riods, which is comparable to the one of the new front 
formation. These remnants show that the observed ref- 
ormation is not an isolated event but part of a quasi- 
periodic reformation cycle. 

The reformation process appears quite coherent along 
the shock surface at the spacecraft separation scale, ~ 
lOAp, thanks to the wavevector alignment with the shock 
normal. Some differences in the pulse shape between CI 
and C2 measurements are nevertheless observed. The 
cause of this inhomogeneity could be the flnite transverse 
correlation length of the upstream ULF wave, which is 
typically of the order of IRe according to Le and Russell 
[ig | or between 8 and 18 Re according to Archer et al. 
Il[. This is however a bit large to account for the inho- 
mogeneity, even when projected on the shock front. It 
is nevertheless possible that the ULF waves are unstable 
to transverse perturbations above some threshold, which 
would reduce their perpendicular correlation length as 
their amplitude grows. Another reason could be fluctua- 
tions or ripples of the shock surface and the inhomogene- 
ity of the ion reflection/ wave amplification process. The 
scale is indeed more consistent with the ripple wavelength 
(a few tens of Ap) estimated by Moullard et al. [25] at the 
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quasi-perpendicular bow shock. In any case this inhomo- 
geneity along the shock front is small since the differences 
between CI and C2 observations are much smaller than 
differences with the other two spacecraft. To a first ap- 
proximation the reformation thus occurs coherently along 
the shock front on scales of at least lOAp. 

V. CONCLUSION 

The global picture is thus in general agreement with 
the scenario of quasi-parallel shock reformation proposed 
by Burgess @: as some of the specularly reflected ions 
escape far upstream they generate waves on the magne- 
tosonic branch which are then convected back into the 
shock and destabilize it. The major difference between 
the Id simulations which inspired this scenario and the 
observations lies in the changes in the wave properties as 
it nears the shock, where wave refraction, an additional 
localized instability mechanism such as the interface in- 
stability and probably nonlinear phenomena participate 
in deflecting, slowing down and amplifying the wave cy- 
cle which ultimately forms a new shock front. The wave 
refraction results in a more coherent reformation along 
the shock surface. On a longer time-scale, this reforma- 
tion process has no reason to stop unless there is a major 
change in the upstream conditions and may be repeated 
quasi-periodically. In this picture the elementary refor- 
mation cycle consists of a shock front and wave locally 
receding at the same speed while the pulse nearest to the 
"old" shock grows to actually become the "new" shock 
front upstream of the previous one, the ramp effectively 
making a jump forward. The old shock slowly vanishes 
and a new pulse grows upstream of the new one, repeat- 
ing the reformation cycle. 

Finally one may list two key elements which conspire 
to create this particular shock structure and dynamics 
lending themselves unusually well to analysis. 

First, there is the shock angle Obu ~ 45°. It means 
that relatively small changes in this angle have impor- 
tant consequence on the kinematics of the reflected ions, 
which either cross the shock after half a gyration and 
thermalise downstream as in quasi-perpendicular shocks 
or escape upstream and feed ULF waves as in quasi- 
parallel shocks. Since on average very few of the reflected 
ions are expected to escape upstream, this angle may be 
the reason why there is such a large gradient in the re- 
flected ion density near the shock front and therefore why 
the wave undergoes such large changes in this region. In 
some sense this leads to a simpler and compacted version 
of lower-angle shocks. For instance the growing pulse 



shares some characteristics with SLAMS [s^] ■ Indeed ev- 
idence suggests that SLAMS grow very rapidly, picking 
up energy from shock-produced energetic ions, and very 
possibly only in close proximity to the shock 23] . In ad- 
dition, they are well-defined entities that propagate from 
upstream into the nominal shock and give rise to the ref- 
ormation of the local shock transition. As such, they 
are manifestations of the dynamic, structured nature of 
coUisionless shocks under oblique geometries. However 
their growth mechanism differs from that of the pulse 
described in this paper, in that they are found in more 
turbulent ULF wave fields and grow off of gradients in dif- 
fuse ion populations rather than the quite mono-energetic 
ions observed here. 

Second, there is the low f3i. A consequence of this is 
that specularly reflected ions form a relatively cold beam 
well distinct from the background ions, helping to shape 
unstable ion distributions. But the low (3 also makes the 
specular reflection itself more delicate by making it more 
difficult for the shock to adjust the electric potential in 
order to reflect the right amount of ions required for the 
dissipation of the incoming solar wind kinetic energy. In 
the zero temperature limit specular reflection even re- 
sults in a " all or nothing" situation which clearly cannot 
be steadily sustained, similarly to the 100% ion reflec- 
tion or transmission occuring in the high Mach number 
perpendicular shock simulations of [26| . The excess of 
reflected ions near the shock front is again a likely possi- 
bility to explain the wave changes in this region, leading 
to the formation of a new shock front. 

These two factors show that the shock is in a situ- 
ation where small changes in Osn or electric potential 
have major consequences, making the shock unstable and 
poorly adaptable to changes in the incoming plasma. As 
it is perturbed major variations in reflected ions strongly 
modify the energy flow and structure of the shock. The 
dynamical modifications of the shock scales and structure 
may then be viewed as a first order attempt by the shock 
to correct an inappropriate reflection of ions while try- 
ing to maintain on average the asymptotic downstream 
state. 



Acknowledgments 

The authors thank M. Fujimoto, S. P. Gary, T. Hor- 
bury and I. Shinohara for fruitful discussions. Parts of 
this work were supported by a grant to Imperial College 
from the UK Science and Technology Facilities Council. 
BL also acknowledges partial support from DOE grant 
DE-FG02-07ER54941. 



[1] Archer, M., T. S. Horbury, E. A. Lucek, C. Mazelle, 
A. Balogh, and I. Dandouras (2005), Size and shape of 
ULF waves in the terrestrial foreshock, J. Geophys. Res., 



110, A05,208, doi:10.1029/2004JA010791. 
[2] Balogh, A., et al. (2001), The Cluster Magnetic Field In- 
vestigation: overview of in-flight performance and initial 



11 



results, Annales Geophysicae, 19, 1207-1217. 

[3] Biskamp, D., and H. Welter (1972), Structure of the 
earth's bow shock., J. Geophys. Res., 77, 6052-6059, doi: 
10.1029/ JA077i031p06052. 

[4] Blanco-Cano, X., N. Omidi, and C. T. Russell (2006), 
Macrostructure of coUisionless bow shocks: 2. ULF waves 
in the foreshock and magnetosheath. Journal of Geo- 
physical Research (Space Physics), 111, 10,205 — h, doi: 
10.1029/2005JA011421. 

[5] Burgess, D. (1989), Cyclic behavior at quasi-parallel col- 
lisionless shocks, Geophys. Res. Lett, 16, 345-348, doi: 
10. 1029/GL016i005p00345. 

[6] Burgess, D. (1995), Foreshock-shock interaction at col- 
lisionless quasi-parallel shocks. Advances in Space Re- 
search, 15, 159-169, doi:10.1016/0273-1177(94)00098-L. 

[7] Burgess, D., et al. (2005), Quasi-parallel Shock Structure 
and Processes, Space Science Reviews, 118, 205-222, doi: 
10.1007/sll214-005-3832-3. 

[8] Eastwood, J. P., A. Balogh, M. W. Dunlop, and T. S. 
Horbury (2002), Cluster observations of fast magne- 
tosonic waves in the terrestrial foreshock, Geophys. Res. 
Lett, 29, 2046, doi:10.1029/2002GL015582. 

[9] Eastwood, J. P., A. Balogh, E. A. Lucek, C. Mazelle, and 
I. Dandouras (2005), Quasi- monochromatic ULF fore- 
shock waves as observed by the four-spacecraft Cluster 
mission: 2. Oblique propagation. Journal of Geophysi- 
cal Research (Space Physics), 110, 11,220, doi:10.1029/ 
2004JA010618. 

[10] Escoubet, C. P., M. Fcliringcr, and M. Goldstein (2001), 
The Cluster mission, Annales Geophysicae, 19, 1197- 
1200. 

[11] Gary, S. P. (1991), Electromagnetic ion/ion instabilities 
and their consequences in space plasmas - A review. Space 
Science Reviews, 56, 373-415, doi:10.1007/BF00196632. 

[12] Gosling, J., M. Thomsen, S. Bame, W. Feldman, 
G. Paschmann, and N. Sckopke (1982), Evidence for 
specularly reflected ions upstream from the quasi-parallel 
shock, Geophys. Res. Lett, 9(12), 1333-1336. 

[13] Gustafsson, G., ct al. (1997), The Electric Field and 
Wave Experiment for the Cluster Mission, Space Science 
Reviews, 79, 137-156. 

[14] Hada, T., C. F. Kennel, and T. Terasawa (1987), Excita- 
tion of compressional waves and the formation of shock- 
lets in the earth's foreshock, J. Geophys. Res., 92, 4423- 
4435, doi:10.1029/JA092iA05p04423. 

[15] Hoppc, M. M., and C. T. Russell (1983), Plasma rest 
frame frequencies and polarizations of low-frequency up- 
stream waves: ISEE 1 and 2 observations, J. Geophys. 
Res., 88, 2021-2028. 

[16] Horbury, T. S., et al. (2001), Cluster magnetic field ob- 
servations of the bowshock: Orientation, motion and 
structure, Annales Geophysicae, 19, 1399-1409. 

[17] Johnstone, A. D., et al. (1997), Peace: a Plasma Elec- 
tron and Current Experiment, Space Science Reviews, 79, 
351-398. 

[18] Krasnosolskikh, V. V., B. Lembege, P. Savoini, and 
V. V. Lobzin (2002), Nonstationarity of strong collision- 
loss quasiporpcndicular shocks: Theory and full particle 
numerical simulations. Physics of Plasmas, 9, 1192-1209. 

[19] Le, G., and C. T. Russell (1990), A study of the coherence 
length of ULF waves in the earth's foreshock, J. Geophys. 
Res., 95, 10,703-10,706, doi:10.1029/JA095iA07pl0703. 

[20] Lcfcbvro, B., S. J. Schwartz, A. F. Fazakcrloy, and 
P. Decreau (2007), Electron dynamics and cross- 



shock potential at the quasi-pcrpcndicular Earth's bow 
shock, J. Geophys. Res., 112, A09,212, doi:10.1029/ 
2007JA012277. 

[21] Lembege, B., et al. (2004), Selected Problems in 
Collisionless-Shock Physics, Space Science Reviews, 110, 
161-226, doi:10.1023/B:SPAC.0000023372.12232.b7. 

[22] Lobzin, V. V., V. V. Krasnoselskikh, J.-M. Bosqucd, 
J.-L. Pincon, S. J. Schwartz, and M. Dunlop (2007), 
Nonstationarity and reformation of liigli-Macli-number 
quasiporpcndicular shocks: Cluster observations, Geo- 
phys. Res. Lett, 34, 5107, doi:10.1029/2006GL029095. 

[23] Lucek, E. A., T. S. Horbury, I. Dandouras, and H. Reme 
(2008), Cluster observations of the Earth's quasi-parallel 
bow shock. Journal of Geophysical Research (Space 
Physics), 113, 7, doi:10.1029/2007JA012756. 

[24] Mazelle, C, et al. (2003), Production of gyrating ions 
from nonlinear wave-particle interaction upstream from 
the Earth's bow shock: A case study from Cluster-CIS, 
Planetary and Space Sciences, 51, 785-795, doi:10.1016/ 
S0032-0633(03)00107-7. 

[25] Moullard, O., D. Burgess, T. S. Horbury, and E. A. Lucek 
(2006), Ripples observed on the surface of the Earth's 
quasi-perpendicular bow shock, J. Geophys. Res., Ill, 
9113, doi:10.1029/2005JA011594. 

[26] Quest, K. B. (1985), Simulations of high-Mach-number 
coUisionless perpendicular shocks in astrophysical plas- 
mas, Physical Review Letters, 54, 1872-1874, doi: 10. 
1103/PhysRevLett.54.1872. 

[27] Quest, K. B. (1988), Theory and simulation of collision- 
less parallel shocks, J. Geophys. Res., 93, 9649-9680, doi: 
10. 1029 / J A093iA09p09649. 

[28] Rome, H., et al. (2001), First multispacccraft ion mea- 
surements in and near the Earth's magnetosphere with 
the identical Cluster ion spectrometry (CIS) experiment, 
Annales Geophysicae, 19, 1303-1354. 

[29] Scholer, M., and D. Burgess (2007), Whistler waves, core 
ion heating, and nonstationarity in oblique coUisionless 
shocks. Physics of Plasmas, 14(7), 072,103, doi:10.1063/ 
1.2748391. 

[30] Scholer, M., M. Fujimoto, and H. Kucharek (1993), Two- 
dimensional simulations of supercritical quasi-parallel 
shocks: upstream waves, downstream waves, and shock 
re-formation, J. Geophys. Res., 98, 18,971, doi: 10. 1029/ 
93JA01647. 

[31] Schwartz, S., M. Thomson, and J. Goshng (1983), Ions 
upstream of the Earth's bow shock: a theoretical compar- 
ison of alternative source populations, J. Geophys. Res., 
88 {A3), 2039-2047. 

[32] Schwartz, S. J. (2000), Shock and discontinuity nor- 
mals, mach numbers, and related parameters, in Analysis 
Methods for Multi- Spacecraft Data, ISSI Scientific Re- 
port, vol. SR-001, edited by G. Paschmann and P. Daly, 
pp. 249-270, Kluwcr. 

[33] Schwartz, S. J. (2006), Shocks: Commonalities in Solar- 
Terrestrial Chains, Space Science Reviews, 124, 333-344, 
doi:10.1007/sll214-006-9093-y. 

[34] Schwartz, S. J., and D. Burgess (1991), Quasi-parallel 
shocks - A patchwork of three-dimensional struc- 
tures, Geophys. Res. Lett, 18, 373-376, doi:10.1029/ 
91GL00138. 

[35] Song, P., and C. T. Russell (1999), Time Series Data 
Analyses in Space Physics, Space Science Reviews, 87, 
387-463, doi: 10.1023/A: 1005035800454. 

[36] Stone, E. C, et al. (1998), The Advanced Composition 



12 



Explorer, Space Science Reviews, 86, 1-22, doi:10.1023/ 
A:1005082526237. 
[37] Thomas, V. A., D. Winskc, and N. Omidi (1990), 
Re-forming supercritical quasi-parallel shocks. I - One- 
and two-dimensional simulations, J. Geophys. Res., 95, 
18,809-18,819. 

[38] Thomsen, M. F., J. T. Gosling, S. J. Bame, T. G. On- 

sagcr, and C. T. Russell (1990), Two-state ion heating at 
quasi-parallel shocks, J. Geophys. Res., 95, 6363-6374, 
doi: 10. 1029/ J A095iA05p06363. 



[39] Thomsen, M. F., J. T. Gosling, T. G. Onsagcr, and C. T. 
Russell (1993), Ion and electron heating at the low-Mach- 
number, quasi-parallel bow shock, J. Geophys. Res., 98, 
3875-3888, doi: 10. 1029/92 J A02560. 

[40] Winske, D., V. A. Thomas, N. Omidi, and K. B. Quest 
(1990), Re-forming supercritical quasi-parallel shocks. II 
- Mechanism for wave generation and front re-formation, 
J. Geophys. Res., 95, 18,821-18,832. 



